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inflation which results from competing boundary condition proposals for the 
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I. INTRODUCTION 

Over the past year new developments concerning the possibility of open inflation P]|| 
in quantum cosmology have led to renewed interest in the debate |2|-f4| about the relative 
merits of competing proposals for boundary conditions on the wave function of the Universe. 
This debate has centred on the differing predictions of the "no boundary" proposal of Hartle 
and Hawking j5|, and the "tunneling" proposals of Vilenkin ||[7] and of Linde [BJ, which are 
themselves distinct. 

Much of the debate between proponents of the competing proposals arises from differ- 
ing predictions obtained using probability amplitudes which are assumed to apply to the 
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respective wave functions. In particular, the nucleation probability for instanton-dominated 
transitions which describe a universe which "tunnels from nothing" is assumed to be 
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where the subscripts (NB), (TL) and (TV) refer to the no-boundary wave function and the 
tunneling wave functions of Linde and Vilenkin respectively. For the solutions in question, 
which correspond to a model in which gravity is coupled to a scalar field, 0, with potential, 
V(<f>) in dimensionless units (see @ below for our conventions), the Euclidean action of the 
instanton is 
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4> being the value of the scalar field at nucleation. 

It has been pointed out by supporters of Vilenkin's proposal P,rC] that arguments which 
use (HD as their starting point are sometimes ill-founded, since the identifications assumed in 
(H|) are not supported by more detailed minisuperspace calculations in models such as that 
of black hole nucleation in the very early universe @,|Kj. In particular, it appears that \I/ TV 
cannot be identified with the probability amplitude (|l|) in this case, thereby side-stepping the 
criticism of ||10| . We have pointed out recently [nj, however, that similar criticisms apply to 
the model involving the nucleation of a closed Friedmann- Robert son- Walker universe with 
an inflationary scalar field - the "work horse model" of numerous discussions in quantum 
cosmology. 
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The conclusion drawn in 



was that the identification of ^ TV , (or indeed of ^ TL ) 



with the probability amplitude ([!]) depends on Planck scale physics through the choice 
of operator ordering in the minisuperspace Wheeler-DeWitt equation. This would not be a 
serious problem if there was a natural choice of operator ordering for which the identification 
([!]) could still be made in the case of ^ TV . However, the choice of ordering which to our 
knowledge is the only one to have been claimed as natural in this minisuperspace model 
I2|-|l4| , is one for which \I/ TV cannot be normalized, the choice of normalization being central 
to the identification ([[]). The grounds for "naturalness" of the relevant factor ordering include 
1-loop unitarity [|TJ]. In fact, the same operator ordering was assumed in calculations which 
claim to lead to the suppression of values of the scalar field with a potential, V(4>), above the 
Planck scale in the relevant probability amplitude [JOT - or even to lead to an enhancement 
of sub-Planckian values of V(4>) in a conformally coupled model [IB] - arguments which 
have usually been assumed 0,[15|,|16|] to favour Vilenkin's tunneling wave function, ^/ Ty . Our 
observation [11], which concerns Planck scale corrections in the scale factor, a, rather than 
in the scalar field <ft (as assumed in ]14-|l6||), casts considerable doubt on such claims. Thus 



the "prediction of inflation" cannot be claimed as a success of Vilenkin's wave function, as 
is often widely assumed. 
The conclusions of ref. 



nj] were based on an analysis of the second order ordinary 



differential equation obtained in the "de Sitter minisuperspace approximation" for arbitrary 
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ordering of the scale factor in the Wheeler-DeWitt equation. Our analysis was based on 
the qualitative behaviour of this equation in various relevant limits. In fact, exact series 
solutions may be developed for the models of interest. It is our intention to present these 
solutions in this paper, and to relate the solutions to the results of [II]. 



II. MINISUPERSPACE MODEL 

The specific model considered here is the 2-dimensional minisuperspace corresponding 
to the classical action for gravity coupled to a scalar field, 
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AtiG = 47rm pl ^ nck , K is the trace of the extrinsic curvature, and the metric is 



where k 2 

assumed to take the closed Friedmann-Robertson- Walker form 



ds 2 = a 2 {-N 2 dt 2 + a 2 (t)dQ 3 2 } , (4) 

where df2 3 2 is a round metric on the 3-sphere, and a 2 = k 2 /(6ti 2 ). 

The Hamiltonian constraint obtained from the (3+l)-decomposition of the field equations 
may be quantized to yield the Wheeler-DeWitt equation 
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where 



U{a,<j>) = l-a 2 V(<f>), 



(6) 



and we have allowed for possible operator-ordering ambiguities through the integer power, 
p, in the first term. The approximation that has been adopted in discussions such as that 
of the prediction of inflation from quantum cosmology |],|7] is to confine the discussion to 
regions in which the the potential V(4>) can be approximated by a cosmological constant, so 
that the <ft dependence in (Bp can be effectively ignored and a standard 1-dimensional WKB 
analysis applied. 

The Wheeler-DeWitt equation for the "de Sitter minisusperspace approximation" has 
been solved exactly in terms of Airy functions by Vilenkin for the special case of operator 
ordering p — —1. We wish to observe here that similar series solutions can be readily 
obtained for the case of arbitrary factor ordering. In particular, with the redefinition 



Z -(P-Wy( z ) 



(7) 



where z = |a 2 , we find that the 1-dimensional Wheeler-DeWitt resulting from (BJ) when the 
^-dependence is neglected is given by 
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z2 Pi + + ( 2Vz3 -z 2 -v 2 )y = o, (8) 

dz z dz v 7 

where u — \-r{p — 1)|. This equation is classified as a type {25} Bocher equation in the ter- 
minology of Moon and Spencer |DJ , although its solution has not been explicitly tabulated. 
Near the regular singular point at z — 0, the solution y(z) approaches that of the modified 
Bessel equation, since for finite V the term involving V in (|]) is sub-dominant. 



III. CONVERGENT SERIES SOLUTIONS 

It is straightforward to apply the Frobenius method to (§) to obtain series solutions, 
although the calculations are rather lengthy. In particular, after a considerable amount of 
algebra we find that the two linearly independent solutions may be constructed in terms of 
the power series y+(z) and y~(z) given by 

y ± {z) = jl + £ [A±z™ - A± +1 z^} } (9) 

where the coefficients of the even and odd powers of z are respectively 

[1/3] 

^ = £( 2 ^f s±% (10) 

i=0 



[(«-l)/3] 

^ + i= E (2^f +1 s3\ (ii) 

j=0 

j ] denotes the integer part, and the sums S± £ are defined iteratively by 

S° ±/ = e° ±i2 „ £ > 0, (12) 

*±4 = E ^S +1 *±j. *>1, <•<[(/ - l)/3], (13) 

J=3i 

*±^= E e ±,^A ^>3, 0<i<[J/3], (14) 

i=3i-l 

(15) 

in terms of the quantities & ± e , defined by 
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j and being two positive integers which differ by an even number. For the special cases 
z/=l,|, 2, the quantities ( |P7D must be replaced by 



- (_i)(*-i)/2 r (| + 1) r (z/ 



2*-ir (| + i) r {v - i) 
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if < 2z/; or k > 2v and 2v — k is odd, 

if /c = 2z/; or k > 2v, j < 2v and 2v — k is even, 

if > 2z/, j > 2v and 2z/ — k is even, 



The linearly independent solutions to @ are then 

*i(*) = y+(*) 



(18) 
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and 
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+ B_y*(z), i/ = 1,|,2,|,.. 



where 
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(with Aq = 1, A± x = 0), and the series y*(z) is defined by 
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The quantities in the summation in the coefficient of z n+u in ( P^D are defined in terms of 
the sequences of all partial sums obtained from partitioning the integer n by 2 or 3. In 
particular, let 

V n ,i = | a} G Z | ai < a{ < . . . < . = n, aj G {2, 3}, (arj +1 - a}) G {2, 3}| (23) 

where X> n> j = dim T 7 ^. The index i = 1,2, . . . ,h n labels the different possible distinct 
ordered sequences of such integers, the number of which is given by 
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[n/3], [n/2] 



E 
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3r+2s=n 
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i.e., the nth element of the Padovan sequence Hl8ipJ| , {1, 0, 1, 1, 1, 2, 2, 3, . . .}, defined recur- 
sively by h = h 2 = 1, hi — 0, h n — h n _ 3 + h n _ 2 for n > 3. The sequences in V n may be 
constructed iteratively by simply appending the integer {n} to each sequence in P n _3 and 
in V n -2, as shown in Table 1. 

For p — 1, (i.e., v — 0), for example, we find solutions given by ([H]) and (|2~0"D with 
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TABLE I. Integer sequences contributing to y* 

71 H n T^n^i't *i> — 1) • • • i 

2 1 {2} 

3 1 {3} 

4 1 {2,4} 

5 2 {2,5}, {3,5} 

6 2 {3,6}, {2,4,6} 

7 3 {2,4,7}, {2,5,7}, {3,5,7} 

8 4 {2,5,8}, {3,5,8}, {3,6,8}, {2,4,6,8} 

9 5 {3, 6, 9}, {2, 4, 6, 9}, {2, 4, 7, 9}, {2, 5, 7, 9}, {3, 5, 7, 9} 

10 7 {2, 4, 7, 10}, {2, 5, 7, 10}, {3, 5, 7, 10}, {2, 5, 8, 10}, {3, 5, 8, 10}, {3, 6, 8, 10}, {2, 4, 6, 8, 10} 
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IV. ASYMPTOTIC SERIES AND THE WKB LIMIT 



Although the series solutions of Sec. Ill converge for all positive values of a, the rate of 
convergence is extremely slow and for many purposes a knowledge of the asymptotic series 
is desirable. For example, boundary conditions are often formulated in terms of the WKB 
solutions, which are effectively the leading order terms in such limits. Since the nature of 
the WKB solutions varies according to the sign of U(a, 0) (Q), instead of using the variable 
z of Sec. Ill, it is convenient to introduce a new variable 



x 



(2Vy 2/3 (1 - 2Vz) 



(27) 



and to make the transformation y 
lently as 



~ x l 2 u. We then see that (Rf) may be written equiva- 
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x H — 



(2V) 4/3 /i 



{2V) 2/3 x-l 



u = 



(28) 



where 



i(p + l)(p_ 3 ). 



(29) 



The essential point to note is of course that to leading order the asymptotic series derived 
from (|28| ) as x —>■ ±oo are no different to those of the special case \i — 0, or p — — 1, 3, for 
which u(x) is an exact Airy function [^]. 

The first limit of interest is x — > oo, which, when physically relevant, corresponds to 
large finite values of a such that a?V <C 1 nonetheless. This requires us to take V -C 1, i.e., 
to restrict the scalar field, 0, to values with a potential much below the Planck scale, which 
is of course physically justified. 

Since \1/ = (|a 2 ) P u(x), we find that the appropriate asymptotic series for \1/ are 
linear combinations of the two modes 
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The leading order term is recognized as the WKB mode in the tunneling region. In practice, 
the requirement that a 2 V 1 means that this expansion is of little relevance beyond the 
leading term: terms of order V in the numerator occur for arbitrarily large odd half-integral 
powers of U in the denominator, and similarly for terms of order V 2 . 

The question as to which linear combination of the modes ( |30D correspond to the re- 
spective solutions \1> 1 = a^'^^Y^a 2 ) and ^ 2 
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expanding them in powers of 2Vz = a V . In particular, from 



of (|5[) can be resolved by 
J)-(f22]) we find Y^z) and 



YJz) are given by their definitions (|19|) and (^) where now 
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(34) 



for large a to compare the leading terms with (|30|) , noting that the overall exponential 
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in (j3~0D may be expanded as exp |4 =p |a 2 (1 + 0(a 2 V)) . In this way we see, 



arly to the observation in fill , that the leading order terms agree if we identify 
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An overall multiplicative factor exp (§^7) must be included in the constant of proportionality 
in fl3U| ) to cancel the overall dependence in the cases that the wave function can be 
normalized. 

In the limit x — > —00, which corresponds to a 2 V 1, or a — > 00, we may find asymptotic 
series which apply beyond leading order. In particular, the linearly independent modes are 
found to be 
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where 
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The phase in (|37]) and (|38|) is obtained by analytic continuation of the corresponding modes 



in (p0|), taking care with branch cuts similarly to the \i = pure Airy function case j20 
To leading order, this of course agrees with the results obtained using the WKB connection 
formulae [11|. The no boundary wave function, ^ NB , of Hartle and Hawking corresponds to 
whereas Vilenkin's tunneling wave function, ^ TV , corresponds to for the 

modes given by ([37]) 



V. DISCUSSION 

The physical interpretation of the series solutions to the de Sitter minisuperspace 
Wheeler-DeWitt equation (^), (0) yield no surprises in comparison with the conclusions 
which have already been reached in fly]]. In particular, for operator orderings with p > 1 
the contribution from Y 2 (^a 2 ) diverges at the origin, causing \1/ to diverge also. (For p < 1 
the divergence is regulated by the prefactor.) Thus any boundary condition proposal which 
requires regularity of the wave function as a — > must reduce to (|35|) with a = if p > 1 
in this particular model; i.e., it must coincide with the proposal of Hartle and Hawking to 
the level of the approximations which have been assumed. Since normalizability of the wave 
function at a = is an essential ingredient of calculations of the probability of inflation 
from quantum cosmology, the question of operator ordering cannot be ignored in such cal- 
culations. If the factor ordering p — 1 is the natural one in the present model, as has been 
claimed [13,14], then there are substantial problems for ^/ TV . While the Hartle-Hawking 



wave function, ^/ NB , has not been found to yield a definitive prediction of inflation for any 
potentials, V(<f>), considered to date pl| , the prognosis for \1/ TV seems to be even more 
doubtful as an identification of type ([I]), (0) cannot be made in the first place if p > 1. 

One might wonder whether these conclusions would be altered by the addition of non- 
inflationary matter degrees of freedom to the model. Although strictly speaking the matter 
degrees of freedom themselves should be quantized, some indication of their bulk effect might 



be obtained by adding a perfect fluid source to the Wheeler-DeWitt equation p2|-Pq] . For a 
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homogeneous isotropic fluid of energy density, p, and pressure obeying an equation of state 
P = wp, we find that the minisuperspace Wheeler-DeWitt equation (|8|) is replaced by 

z *^y + z ^- + hVz 3 + (3z 3 ^/ 2 - kz 2 -u 2 )y = 0, (41) 
dz 2 dz v ' 

where the spatial curvature, k, is identified explicitly for the purpose of the discussion that 
follows, z = \o? and v = |(p — 1) as before. The normalization factor in (]3|) is now given in 
general by a 2 = k 2 / (3V), V being the 3- volume of the unit spatial hypersurface (with suitable 
topological identifications for k = 0,-1), and (3 = 2n 2 a 2 a^ 1+w ^ /(3p Q ) for a matter source 
with energy density normalized to p = p at a = a Q . For large values of a, the scalar potential 
term dominates irrespective of the other matter degrees of freedom or the spatial curvature, 
and the solutions approximate to a superposition of modes a^~ p ^ 2 J( p -i)/§{\a 3 '\/V) and 



-i)/§(\a?\/V) as was discussed in [[□] 



Now let us consider the properties of the solutions as a — > 0, which is important for 
normalization of the wave function. The case of stiff matter (w — 1) in general yields a 
singular wave function as a — > for both modes, and hence any choice of boundary condition, 
the singularity being of infinitely oscillatory type if (3 > v 2 = -^(p — l) 2 , 

^ a ~ Q a^' 2 \C 1 cos lnd^-z/ 2 a 2 ) + C 2 sin ln^^-i/ 2 a 2 ) j, (42) 

where C 1 2 are constants. This is a pathological type of behaviour which does not distinguish 
between different boundary conditions for the wave function, and is presumably not of 
physical interest. Apart from this case, however, one may readily see that the conclusions 



of [|TTJ are not altered in general by the presence of additional matter, since the indicial 
equation for series solutions about a = 0, which primarily determines the properties of the 
solutions near a = 0, is unaffected for w < 1. 

Consider the case of radiation (w = |), for example, in which case ([|l|) is still of Bocher 
type {25} [0. As a — > the leading order behaviour is given by that of the exact solutions 
to ( f4"l"D with V = 0, which are given by P5J 

* = a -^ +1 ) [C x M m ^ 2 ) + C 2 W mu (a 2 )} , 

= exp[-ia 2 ] {C.M^ + | - f , 1 + 2u, a 2 ) + C 2 \j(v + \ - f , 1 + 2u, a 2 ) } , (43) 

in terms of Whittaker functions and Kummer functions |27]] respectively, for 
k = +1. Since exp[-ia 2 ]M(z/ + |, 1 + 2v, a 2 ) = I> + 1] I u (\a 2 ) and 



exp[— ^a 2 ]u(^u + |, 1 + 2u, a 2 ^j = ■n~ x l 2 ar 2v K u (^a 2 ) [p7|1 , we of course retrieve the limit- 
ing V = solutions discussed in previous sections when j3 = 0. The solutions for k = — 1 
and V = are likewise given by 

$ = l Cl [e ia2 ' 2 M(v + | - if, 1 + 2z/, -m 2 ) + e - w2 / 2 M(z/ + ± + if, 1 + 2u, ia 2 ) } 

+ \C 2 {e-^ +ia2 / 2 u(u + | - if , 1 + 2z/, -ia 2 ) + e l7Tl/ ~ ia2/2 u(u + | + if , 1 + 2z/, ia 2 ) } 

(44) 
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which reduces to a superposition of ordinary Bessel function modes 



= (i-p)/a [b^q 2 ) + B 2 Y u {\a 2 )} , (45) 
when /3 = 0, where relative to (gD 5 X = 2 2 *T[z/ + 1]^ and 5 2 = — |VttC 2 . In the spatially 



flat case, k — 0, the solution to (¥1) with V = is given in terms of ordinary Bessel functions 
modes for all values of (3 p4| , p6 | 



* = a (1 ~ p)/2 \c x J 2v {Jpa) + C 2 Y 2 ,(V^«)} • (46) 



In each of the cases (j43|) , (|44T) and (f46|), the mode with coefficient C x is finite as a — > 0, 
whereas the mode with coefficient C 2 is proportional to a~ 4u = a l ~ p as a — > 0, which causes 
\I/ to diverge for p > 1 unless C 2 = 0, in direct analogy to the radiation-free solutions 
(P = 0). Unfortunately, the requirement that Vilenkin's wave function, \I/ TV , is outgoing 
means that the coefficients of both modes must be non-zero, and C 2 = is not possible for 
\1/ TV . In particular, in the model here Vilenkin's boundary condition requires that 

?' <9\I> 

Re > °- < 47 > 

in the Lorentzian region, which includes all values a > if k = 0, —1. In the k = case 
(|46|), for example, the linear combination of modes which satisfies this condition is 

^ TV oc a^' 2 | J 2u (J(3a) - iY 2v {^a)\ , (48) 

namely a Hankel function of the second kind, for which Re ^^ 8 ^ v = 2[ira( J2J 2 + ^u 2 )]^ 1 - 
Thus the problems faced by Vilenkin's wave function, especially with regard to the prediction 
of inflation |TTJ, appear to be generic. 



One must also conclude that great care must be exercised when making arguments based 
on the identification of probabilities according to ([!]), (@). In general, the contribution of 
the prefactor and choices which depend on Planck scale physics cannot be ignored when one 
is considering the problem of the nucleation of the universe. 
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